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Abstract 

Glauber gluons in Drell-Yan processes are soft gluons with the transverse momenta much 
larger than their momentum components along the directions of initial hadrons. Their existence 
has been a serious challenge in proving the factorization of Drell-Yan processes. The recently 
proposed soft collinear effect theory of QCD can provide a transparent way to show factorizations 
for a class of processes, but it does not address the effect of glauber gluons. In this letter we 
first confirm the existence of glauber gluons through an example. We then add glauber gluons 
into the effective theory and study their interaction with other particles. In the framework of 
the effective theory with glauber gluons we are able to show that the effects of glauber gluons in 
Drell-Yan processes are canceled and the factorization holds in the existence of glauber gluons. 
Our work completes the proof or argument of factorization of Drell-Yan process in the framework 
of the soft collinear effective theory. 



QCD factorization theorems enable us to make theoretical predictions from QCD for hard- 
scattering processes, where one can separate long-distance effects from short-distance effects. The 
short-distance effects depend on the details of a process and can be studied with perturbative 
QCD, while long-distance effects can be characterized by matrix elements of QCD operators, which 
depend only on the structure of hadrons and infrared properties of QCD[lJ. In general, it is 
difficult to prove a factorization theorem, e.g., much effort has been spent to prove factorization for 
Drell-Yan process [21 [3]. A main obstacle has been the exchange of glauber gluons in initial state 
interactions [U [3]. The glauber gluons are soft with the transverse momenta much larger than their 
momentum components along the directions of initial hadrons. 

Recently, the soft collinear effective theory(SCET) has been proposed to study interactions 
among collinear- and various soft particles of QCD [5j El [7]. The interactions are at long-distance and 
produce collinear- and I.R. divergences in perturbative calculations. Assuming that SCET includes 
all those particles of QCD, then SCET will re-produce these collinear- and I.R. singularities. It 
is interesting to note that with SCET the proofs or arguments of factorization theorems can be 
simplified or made in a transparent way. In p] several examples including Drell-Yan processes 
are given. However, the glauber gluons are not included in the proposed SCET[5l [6]. In this 
letter we make an attempt to include the glauber gluons into SCET and show that their effect in 
the factorization of Drell-Yan process is canceled [6]. Therefore our work completes the proof or 
argument of factorization of Drell-Yan processes in the framework of SCET. 

Before our study of SCET by adding glauber gluons and the factorization of Drell-Yan processes, 
we show the existence of glauber gluons in scattering amplitudes. For this we consider the process 
<l(pi) + <l(p2) ~~ ^ 9(k g ) + 7* (q), which contributes to the differential cross section of Drell-Yan 
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Figure 1: A diagram which contributes to the process q(pi) + q{p2) g(k g ) + 7*(<?) at one loop 
level. The broken line is a cut. The black circle denotes the insertion of the electromagnetic current 
operator. 

processes. At one loop level, the amplitude of the process receives a contribution from Fig. la. We 
will use the light-cone coordinate system, in which a vector is expressed as = (a + ,a~ ,a±) = 
((a + a 3 )/\/2, (a — a 3 ) / y/2, a 1 , a 2 ) and a\ = (a 1 ) 2 + (a 2 ) 2 . We introduce two light cone vectors 
n M = (1, 0, 0, 0) and = (0, 1, 0, 0). The momenta are given as: 

tf = (pt, 0,0,0), p£ = (0,^,0,0), k% = (k+,k~,k g± ). (l) 

We consider the contribution in the momentum configuration where the momentum k g scales as 
kg 1 ~ Q(A 2 ,1,A, A) with Q as a large scale and A <C 1. In this configuration the outgoing gluon 
is collinear to the initial q. For convenience we will take Q = 1 in the following. The leading 
contribution from Fig.l in the limit A — > can come from different momentum regions of the 
loop momentum k with k 2 ~ A 2 or with k 2 at higher order of A. These possible regions are: The 
collinear regions with k at the order of (1,A 2 ,A,A) or (A 2 ,1,A,A), the soft region with k at the 
order of (A, A, A, A), the ultra-soft region with k at the order of (A 2 , A 2 , A 2 , A 2 ) and the glauber 
region with k at the order of (A 2 , A 2 , A, A). Gluons with momenta in glauber regions are called as 
glauber gluons. To show that the glauber gluon gives nonzero contribution to the amplitude, it is 
enough to consider the following loop integral involved in the contribution of Fig. la: 

T- [ ^ 1 1 1 1 (2) 

J (2vr) 4 k 2 + ie {pi - k) 2 + ie (k g - k) 2 + ie (k - k g + p 2 ) 2 + is' {) 

One can perform a power counting in A for the integral in differen momentum regions to determine 
the relative importance of contributions in deferent regions and leading contributions. One finds 
that the contributions from the collinear region with k ~ (A 2 , 1, A, A), the ultra-soft region and the 
glauber region are at leading order of A which is 0(A -2 ). The contributions from other regions are 
at higher order of A. It is complicated to have a clear separation of leading contributions from the 
three regions, because they are all overlapped. 

For our purpose we can consider the absorptive part of Fig. la instead of the total contribution. 
We consider the cut diagram of Figl.a, which is given as Fig. lb. Instead of the integral I we need 
to consider the absorptive- or imaginary part of /. The imaginary part is given by: 

1 r d A k 1 1 

From these on-shell conditions and with ~ (A, A) one always has k + ~ A 2 and k~ ~ A 2 . This 
indicates that the absorptive part of / at the leading order of A is determined only by the glauber 
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region. Performing the power counting for the glauber region one obtains the leading contribution 
to Iml as: 
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with s = 2pfp2- In the above we have used dimensional regularization with e = 4 — d to regularize 
divergences. From the result the integral is divergent. The divergent contribution comes partly 
form the region of k± ~ and partly from the region of k± ~ k g ±. The singular contribution from 
the region of k± ~ k g ± is because that there are pinched poles in the complex A; + -plan for the 
integral /. The singularities generated by these pinched poles have brought up a difficulty to prove 
the factorization in Drell-Yan processes. 

From the above discussion we have Iml = rm/ G (I + 0(A)). To verify this one can work out 
the exact result of the integral /. The exact result can be found in [8J. Through an analytical 
continuation and an expansion of the exact result we have from [8]: 
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Comparing the result in Eq.(4,5) one verifies Iml = LxlZg(1 + C(A))- One can conclude that the 
exchange of glauber gluons does give nonzero contributions to scattering amplitudes at leading order 
of A and hence possibly to differential cross sections. The possible contributions due to glauber 
gluons are at the same order of contributions due to collinear- and ultra-soft gluons. Therefore, one 
should include glauber gluons into SCET as an effective theory of QCD to describe interactions 
between collinear partons. Although the contributions of glauber gluons in Drell-Yan processes at 
twist- 2 level are canceled, as shown in [21 [3] and in the current work, but the existence of glauber 
gluons has observable effects at twist-3 level. Recent studies of single transverse spin asymmetries 
in Drell-Yan- and DIS processes have been shown that these asymmetries are generated through 
exchanges of glauber gluons [9]. Hence it is also important to have SCET containing glauber gluons 
for a possible study of factorizations of single transverse spin asymmetries. 

To derive SCET containing glauber gluons, we first consider SCET of one collinear system, later 
we will add another collinear system for our purpose. We will use the space-time representation 
of SCET as used in [7]. We choose the collinear system moving along the n-direction with the 
momentum p at the order 

^~Q(I,A 2 ,A,A), (6) 

where Q is a larger scale while A is small with A <C 1. Again we will take Q = 1 in the following. 
Interacting particles in such a collinear system can be classified with their momenta at different 
orders of A. They are collinear particles with momenta at the order p^ ~ (1,A 2 ,A,A) , ultra-soft 
particles with momenta at the order p^ ~ (A 2 , A 2 , A 2 , A 2 ), and glauber gluons with momenta at 
the order p^ ~ (A 2 , A 2 , A, A). The later have not been incorporated in SCET. Soft particles with 
momenta at the order of p^ ~ (A, A, A, A) will not interact with the above particles because of 
momentum conservation. However, their effect will appear in matching operators of full QCD into 
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those of SCET. It has been shown that their effect is canceled in Drell-Yan process. We will not 
consider soft particles here. 

For the collinear system we decompose the gluon field and quark field of full QCD as: 

A» = A£ + A% + A» S , ^ = ^ + v + q. (7) 

In the above An and £ n are the collinear fields for collinear gluons and collinear quarks, respectively. 
77 is a small component for collinear quarks which will be integrated out and expressed in terms of 
other fields. A us and q are the ultra-soft fields for ultra-soft particles. A^ is the field for glauber 
gluons. The space-time dependence of these fields are characterized by the small scale A in the 
following patten. For a collinear field 4> n which can be the collinear gluon field or the collinear 
quark field we have 

~ 0(1)^0*0, d-<j> n (x) ~ 0{\ 2 )<j) n {x), d^n{x) ~ 0{\)<j) n {x). (8) 
For a generic ultra-soft field 4> us and the glauber gluon field we have 
d^ us (x) ~ O(A 2 )0 us (x), 

d + A»(x) ~ 0(X 2 )A^x), d-A»(x)~0(\ 2 )A»(x), P ± A${x) ~ 0{\)A^x). (9) 

By inspecting propagators of particles in different momentum modes we can derive the following 
power-scaling: 

At ~ 0(1), Ar~0(A 2 ), ^ ± ~0(A), e-O(A), 

A» us ~ 0(A 2 ), ^~0(A 2 ), g~0(A 3 ). (10) 

Using the power-scaling of fields and the scaling of the space-time dependence of the fields one 
can expand the Lagrangian of full QCD in A with the above fields. The expansion is straightforward. 
To express the leading order results we introduce the modified derivative or projection. The 
derivative only acts on collinear fields, it gives zero if the derivative acts on ultra-soft- or glauber 
gluon fields, i.e., 

dt<j>n{x) =d + 4> n (x), d£(f> us {x)=0, 8±A% = 0. (11) 

The derivative acts as an usual derivative on all fields while <9^j_ acts only on collinear- or 
glauber gluon fields, but not on ultra-soft fields. The leading terms of our effective theory can be 
written as: 

(£, q, A£,A%, A^ s ) = £ g (A%) + C us (q, A&) + C n (t n , A£, n-A g ,n- A us ) + • • • , (12) 

where • • • stand for terms whose contributions to the action S = J d A xC^ are suppressed by A. 
In the expansion of the action S one should assign the space-time integration cftx with a correct 
scaling in A by noticing that the integration should not change the momentum orders of fields in 
different momentum modes. 

At leading order glauber gluons can not interact with ultra-soft particles and they also can not 
interact with themselves. Hence C g (A g ) is only the kinetic term of the glauber gluons: 

C g = -\Ti[d^ g - dlA»][d^A gv - d ±u A m }. (13) 
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It should be noted that only the transverse derivative appears. This fact results in that the 
propagator of glauber gluons with the momentum k will be proportional to l/k\, as already seen 
in Eq.(4). C us (q, A^ s ) is the part for ultra-soft particles: 



C us (q,A^ s ) = qi 7 • D us q - ^Tr {[£>&, D» us }} 2 , 



(14) 



where the covariant derivative is defined only with the ultra-soft gluon field, i.e., D£ s = + igA£ s . 
Cn(^.n, An, A g , A us ) contains collinear particles and their interactions with ultra-soft- and glauber 
gluons. To express it we introduce the covariant derivative with the modified derivative: 

D^ = d^ + igA£ + ign^n ■ (A g + A us ). 

With the covariant derivative the part L n can be expressed as 

£n(£n, A£,n ■ A g ,n- A us ) = l fl in ■ ■ D n £ n - £^7 • D c± -^—^n ■ 7*7 • D c± £ n 



(15) 



~2iD2 



-^Tr{[^,^]} 2 , 



(16) 



where the collinear covariant derivative D c is defined as = <9 M + igA£. The glauber gluon field 
appears in the covariant derivative in the same way as the ultra-soft gluon field. However, the 
interaction of collinear gluons with ultra-soft gluons are different than that with glauber gluons 
because the operator d n ± does not act on the ultra-soft field, but on the glauber field, i.e., d£ ± A g 7^ 
0. 

The effective theory at leading order A is invariant under various gauge transformations, which 
are collinear-, ultra-soft- and glauber gauge transformation. The space-time dependence of these 
transformations is characterized in the same patten as that of the corresponding fields in Eq. (8,9). 
These transformations are: 



Collinear : A^ s 



av 

Ultra - soft : A% 
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Glauber : A» 
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U c A*Ut - -U c \d» + ig(A» s + A*), U\ 



9 



UusAfJJ^s, A^t — > XJ us A!iXJ us , £ n —7- U us ^ n , q —> U us q, 



g ' ^■US^'-g^uS 
^ 1 ^ US "> 

Sn U g £n, q 

I 



U us A!^ s U U s 
9 



Ap. 



ai -> UgAm- -U g er + igiA^ + AfiM 
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(17) 



Under these transformations the effective Lagrangian is invariant at the leading order of A, or 
the change of the effective Lagrangian is suppressed by A. One can re-express these transforma- 
tions with the modified derivative and neglect higher order effects introduced by the fields in the 
transformation so that the effective Lagrangian is exactly invariant. In the above glauber gauge 
transformation the glauber field is not changed. However, there is another type of transformations 
related to the glauber transformation defined by U g j_(x) = U g (x + n+x~n), under which the glauber 
field is changed: 



At* 



U 9 ±A%UI ± 



l -u, 



9^ 



AV — V A 1 * 



U g ±A£ul ± + U g± [A*,, Ul ± ] , & Ugj.Z*. 



(18) 
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Under the above transformation the Lagrangian is invariant. One needs to fix gauges for different 
gauge fields and introduces corresponding ghost fields. We note that the gauge fixing for ultra-soft- 
and glauber gluons can be fixed separately. We will take Feynman gauge for ultra-soft- and glauber 
gluons. For collinear gluons we take the covariant gauge for collinear gluons as in [5], where the 
covariant derivative should include the glauber gluon field as in Eq.(15). 

To study the effects of glauber- and ultra-soft gluons in Drell-Yan process one needs to include 
another collinear fields in the n-direction. This can be done by adding the collinear gluon field 
and the collinear quark field £ n in the decomposition in Eq.(7) correspondingly and by doing the 
expansion in A. At leading order A the collinear fields in different directions can no interact with 
each other. Momentum conservation insures that the collinear fields in different directions can not 
be mixed with each other. Then, the soft collinear effective theory for two collinear systems moving 
in the n-direction and in n direction respectively can be easily written down: 



£(nn) 



+ £-us{Qi ^4«s) "I" £-n{£,m A^, n ■ A g , n ■ A us ) + £ n (£ n , A£, n ■ Ag, n ■ A us ) + • • • . (19) 



From the effective theory, the interaction between the collinear fields in different directions is only 
through exchanging glauber- and ultra-soft gluons. 

Without the presence of glauber gluons the decoupling of ultra-soft gluons can be simply made 
through redefinitions of collinear fields, as shown in [6]. In the presence of glauber gluons in 
SCET, one can also redefine fields so that the redefined fields do not interacting with glauber gluons 
and ultra-soft gluons. We introduce the following SU (3)-matrices: 

Zfi(x) = Pexpj— ig J dfifi ■ Agifin + x)\ , Z n (x) = Pexp j— ig J dpn ■ A g (fin + x)| , 
Y n (x) = Pexpl^-igJ d(3n ■ A us {j3n + x)\ , Y n {x) = P exp i-ig J dj3n ■ A us (f3n + x)| . 



(20) 



With these matrices we redefine all collinear fields as: 



4^ 



Aft 



Zfi 



1 n-n-n 1 n 
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dji + ig s n^A us , Z\ , £ n = Z n Y n xn, 



Yn-A^Xn Z\ ^Z n + igsn^A^s, Z n , ^ n — Z n Y n Xn- 



(21) 



With the re-defined fields the Lagrangian containing the interactions with glauber gluons at the 
leading order of A becomes 

£ n (Z n , A£, n-A g ,n- A us ) = C n ( X n, A£, 0, 0), C n (^ n , A£, n-A g ,n- A us ) = C n ( Xn , A£, 0, 0). (22) 

This result tells that the redefined collinear fields will not interact with glauber gluons and ultra- 
soft gluons. It should be noted that there is no difference between S- matrix elements calculated 
with the original Lagrangian and those calculated with the Lagrangian of redefined fields. For our 
redefinition of fields this can also be realized with path-integral of SCET by noting the fact that 
the Jacobian's of changing field variables are unit. 

Now we are in position to discuss the factorization for Drell-Yan process. The process is given 

as 

K + h n -> £ + + r + X, (23) 
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where one hadron h n moves in the n-direction while another hn moves in the n-direction. At leading 
order of QED the differential cross-section is determined by the hadronic tensor: 

W^(q)= f d A xe-^ x (h n h n \J^(x)r(0)\h n h n ), J» = ^^. (24) 



To show the factorization one needs first to match the operator J^(x)J v {Q) to the operators of the 
soft collinear effective theory at the leading order in A. In the matching, the interactions between 
collinear particles which are not presented in SCET are assembled by different Wilson lines: 

W n (x) = Pexp |— ig J dan ■ A n (an + x)| , W n (x) = Pexp |— ig J dan ■ A n (an + x)| . 

' ' ~* (25) 

At leading order of A the operators consists of collinear fields only and their relative space-time 
dependence can only be of x~n for <fi n fields and of x + n of <p n . Taking the operators consisting 
only of quark fields in SCET as example, the matching is: 

W^ u (q) = J dkfdk^C^iq, kffi, k^n) J dx~dy + exp (-ix'kt - iy + k 2 ^j 

■(h n h n \ [^nWn)(x-n)T^(wkn)(0)} [^nW n )(y + n)T^(W^ n )(0)] \h n h n ) + • ■ (26) 

where • • • stand for other possible operators at leading order and terms at higher orders in A which 
will be neglected. The leading order is at A 4 . For simplicity we will discuss the factorization for 
the above term in detail. Those terms containing other possible operators can be handled in the 
same way. In the above we have already used the properties 

n ■ 7£n = l~£n = 0, n • 7^ n = 7 + £„ = (27) 

to decouple the Dirac indices. Hence the T^-matrices are given by 

(r«, if } , rf ) = (n ■ 7 , n ■ 775, n ■ 77L ) ■ (28) 

Y^n is defined by replacing n with fi in if. 

The functions C can be calculated with perturbative theory. If one uses perturbative theory 
to calculate the hadronic tensor by replacing the hadrons with parton states, one will have I.R.- 
and collinear singularities. These singularities are reproduced by the matrix elements of SCET 
operators through construction of SCET, e.g., those in the right hand side of the above equation. 
Hence the perturbative functions are free from I.R.- and collinear divergence. At this stage the 
long-distance- and short-distance effects are factorized, but the factorization of Drell-Yan process 
is not complete. The factorization is completely achieved if one can show that the matrix element 
of the two-hadron state in Eq.(26) can be factorized into a product of two matrix elements of 
one-hadron state. In Eq.(26) the operator (£nWn)(x~ n)lf (W / s?n)(0) can emit ultra-soft- and 
glauber gluons to interact with the state h n and the emitted gluons can also be absorbed by the 
operator (^ n W n )(y + n)Tn \w^ n )(0) . This prevents us at first look from factorizing the matrix 
element of two-hadron state in Eq.(26) as a product of two matrix elements of one-hadron state. 
For the factorization one needs to show that the total effects of gluon-emissions are canceled. 
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Before discussing the factorization, we discuss a general feature of the A-expansion in the space- 
time representation. We consider a product of collinear fields (f>n, denoted as &n, multiplied with 
any ultra-soft field 4> us . The integral of the product has the following property: 



j dft<$> n {ftn + x)<j) us (f3n + x) = J df3<S> n (ftn + x) <j) us (x) {l + 0{\ 2 )) . 



(29) 



This equation can be derived by expansion the ultra-soft field as (p us (ftn+x) = (j) us {x)+ ftd + <p us {x) + 
■ ■ -. The derivative d + 4> us (x) is at the order of A 2 , while the integration variable ft is related to the 
collinear direction and is taken as at the order of A . Therefore we have the above equation. Now 
we use the redefined collinear fields in Eq.(21) to express the matrix elements in Eq.(26). Under 
the re-definition the Wilson W n and W n become by using the property 



W n (x) = V n (x)W n (x)V r l(x){l + 0(X 2 )}, V n (x) = Z n (x)Y n (x), 
W n (x) = Vn(x)Wn(x)V r l(x) jl + 0(A 2 )} , V n (x) = Z n (x)Y n (x) 



(30) 



where W n and W% are defined by replacing A n and A n in Eq.(25) with A n and A n , respectively. 
The hadronic matrix element with the redefined fields become: 

(hnK\ \{tnWn){x-n)Tf(wU n m] WnW n ){y+ n)T$ (W^ n )(0)] \h n h n ) 



{hnh n \(xnWn)i(x-n)T { i\wlxn) j (0) M(x-n)Vn(0)}.. \v^(y + n)V n (0) 
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ki 



{XnW n ) k {y + n)Y^\wlxn)m\Khn) [l + 0(X 2 ) 



(31) 



where ijlk are color indices. With the definition of Vn(x~n) and V n (y + n) one knows that the 
^"-dependence in Vn{x~n) and the y + -dependence in V n (y + n) are characterized by A 2 . Since the 
space-time separation x~ and y + will be integrated over in the hadronic tensor, hence we can have 
the following approximation under the integrations: 



vi(x-n)V n (0)].. fcV n)V n (0)] = 8 i3 5 lk + 0(X 2 ). 

J tj l J kl 



(32) 



Using the fact that redefined collinear fields do not interact with glauber- and ultra-soft gluons and 
also that the collinear fields in different directions do not interact, we can re-write Eq.(31) under 
the integrations in Eq.(26) with Eq.(32) as: 

(hnK\ [(ZnWn)(x-n)rg>(Wkn)(0)} [{l n W n ){y + n)Y^ (W^ n )(0)] \KK) 
= [{h n \(xnWn)(x~n)Tf(wlxn)(0)\hn)] [(K(XnW n )(y + n)T^\wlxn)(0)\hn) 
•{l + 0(A 2 )}. (33) 

Substituting this result and expressing the redefined fields with the original fields, we complete the 
factorization for Eq. (26): 



W^{q) = j dk+dk^C^(q,k+n,k^n) ■ j dx~ 



(hn I (& Wn) {x-n)Tf iWUn) (0) | K) 



d, 



J dy+e- i y +k ^h n \(CnW n )(y + n)r^\W^n)mhn) 



1 + 0(A 2 



+ 



(34) 
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In the matching in Eq.(26) there are contributions from color-octet operators at leading order like 

(35) 

represented by ■ ■ ■. With the above procedure and the fact that the hadrons are color- less, one can 
show that these color-octet operators do not contribute at leading order of A. 

In the above we have shown that the effects of glauber gluons are completely canceled at the 
leading order of A in the same way of the cancelation of ultra-soft gluons. The above result is 
invariant at leading order of A under various gauge transformation. We note that the Wilson line 
W n and Wfi under the transformations given in Eq. (17) are transformed as: 

Collinear :W n -> U n W n (l + 0{\ 2 )) , W n -> U n W n (l + 0{\ 2 )) . 

Ultra - soft : W n -> U us W n Ul (l + C(A 2 )) , W n -> U us W n Ul (l + 0(A 2 )) . 

Glauber :W n -> \J g W n Jj\ (l + C»(A 2 )) , W n -> UgWnU* (l + C(A 2 )) . (36) 

At the leading order of A there are other three operators appearing in the matching in Eq.(26). 
They are: 

t,nC,n^n > ?n?n<J„ ^ n > U n U n i \ 6 1 ) 

where we have suppressed the color structure, spin structures, space-time dependence and the 
Wilson line of collinear gluon in different directions. The collinear field strength tensors are defined 
with corresponding collinear gluon fields, the indices fj,, u, // and v' are transverse. The factorization 
or the cancelation of effects through exchanges of glauber- and ultra-soft gluons can be shown in a 
similar way as in the above. With the factorization one can define for each hadron various parton 
distributions and write the final result for the hadronic tensor in a compact form. The detailed 
results can be found somewhere in the literature. 

To summarize: In this letter we have studied the effects of glauber gluons in the framework of 
SCET, which had not been considered in the originally proposed SCET. We have confirmed the 
existence of glauber gluons through an example. By adding the glauber gluons into SCET we have 
found that the glauber gluons only interact with collinear particles. Furthermore, we have shown 
in the framework of SCET with glauber gluons that the effects of glauber gluons are canceled in 
Drell-Yan process, i.e., the factorization in the present of glauber gluons still holds. Therefore our 
work completes the proof or argument of factorization of Drell-Yan process in the framework of 
SCET. Our result is in the agreement with the traditional proof of Drell-Yan process where the 
existence of glauber gluon has been a challenging difficulty. However the proof here is given in a 
more transparent way. 
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